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Hydrogen-atom spectrum under a minimal-length hypothesis 
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The energy spectrum of the Coulomb potential with minimal length commutation relations 
[Xi,Pj] — ih{6ij(l + f3P 2 ) + /3'PiPj} is determined both numerically and perturbatively for ar- 
bitrary values of /3' /{3 and angular momenta I. The constraint on the minimal length scale from 
precision hydrogen spectroscopy data is of order of a few GeV -1 , weaker than previously claimed. 
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Quantum gravity incorporates Newton's constant as 
a dimensional parameter that could manifest itself as a 
minimal length in the system. Recent string theoretic 
considerations suggest that this length scale might imply 
an ultraviolet-infrared (UV-IR) correspondence, contrary 
to the normal perceptions on momentum and spatial sep- 
arations. Large momenta are now directly tied to large 
spatial dimensions, which then implies the existence of 
a minimal length. Earlier studies have focused upon its 
amelioration of ultraviolet divergences 0, but did not 
take into full account the UV-IR correspondence. 

There are various ways of implementing such an idea, 
but the simplest is to suppose that coordinates no longer 
commute in D-dimensional space. This, in turn, leads to 
a deformation of the canonical commutation relations. In 
our previous works, we adopted the equivalent hypothesis 
that the fundamental commutation relations between po- 
sition and momentum are no longer constant multiples of 
the identity. In this paper, we report on constraints on 
the minimal length hypothesis from precision measure- 
ments on hydrogenic atoms. This system has a potential 
that is singular at the origin, and is therefore particu- 
larly sensitive to whether there is a fundamental minimal 
length. Considerations based upon higher-dimensional 
theories suggest that such lengths may be large . 

To set the context, we note that if in one dimension 
we have 



[X, P] = ih(l + pP 2 



where [3 is a small parameter, then the resulting uncer- 
tainty relation (AX)(AP) > ih{l + /3(AP) 2 } exhibits a 
form of the UV-IR correspondence, and gives as minimal 
length AX > hyf]3 0. 

We had examined the harmonic oscillator system un- 
der this hypothesis in Q, but no real constraint can be 
obtained on the minimal length, presumably because of 
the softness of the potential at the origin. An interest- 
ing approach is to take the classical limit fi — > of the 
commutation relations; it yields an unbelievably strong 
bound, but its robustness might be questioned [5|- 

We will work in arbitrary D > 1 dimensions, where l[T|l 
takes the tensorial form 



[X h Pj] = i7L{<%(! + (3P 2 ) + 0'PiPj}, 



(2) 



which, assuming that the momenta commute [Pi, Pj] = 0, 
leads via the Jacobi identity to the nontrivial position 



commutation relations 

(2/3 - p) + (2/3 + (3')pP 2 



[X{ , Xj 



ih ■ 



(1 + /3P 2 ) 



(PiXj - PjXi). 



The position and momentum operators can be repre- 
sented by 

Xi = (l + (3p 2 ) Xl + (3'piPjXj , Pi = p h (3) 

where the operators Xi and pj satisfy the canonical com- 
mutation relations [xi,pj] = ihSij. The simplest repre- 
sentation is momentum diagonal, 



r - fi l 

Xi in , 
9Pi 



Pi = Pi- 



(4) 



In this representation the eigenvalue equation for the 
distance squared operator R 2 = A^A^ can be solved ex- 
actly. With 



09 + 13') P 2 



1 



09 + f3')p 2 



(5) 



the eigenvalues r 2 £ = K 2 (f3 + fi')p 2 n n and eigenfunctions 
Rni are given by (see 3 for details) 



p 2 ne = (2n + a + b+ l) 2 - (1 - V ) 2 { L 2 



(D-lf 



(1) Rra(z)<x (/3 + /3') D/4 



1 — z \ ' I 1 + z x 



2 J V 2 
where P^ L,b \z) are the Jacobi polynomials and 



4 



V 



[! + (£>- 1)# 



l + QP-l)r? , 
A = h a. 



7? 2 L 2 , 



(6) 



Having diagonalized R , one can express the action 
of the operator on any function of definite an- 

gular momentum ^(z) — J2^Lo fnRni( z )- In particu- 
lar, the Schrodinger equation for the Coulomb problem, 

can be rewritten in the 

variable z as 



5> 



l + z 



Pnt 



R nl {z) = 0, (7) 



2 



£ = Axmin/ao being the ratio of the minimal length 
Ax m i n = fiy/P + ft 1 to the Bohr radius a = 1/km, and 
e = £ 2 (E/E ) the energy in units of the usual ground- 
state energy Eq = —Xjla^m times £ 2 . 

Using the recursion relation for Jacobi polynomials as 
well as the orthogonality of the distance cigcnfunctions 
R n e, the Schrodinger equation is equivalent to a recursion 
relation for the expansion coefficients, 

lOn-1 = 0, (8) 

with /_! = 0, s n = 1 - e + 2£/p nt , 



n y n> n {2n + a + b){2n + a + b + 2Y 
and 

2 

fl " ~ ~ (2n + a + b + 2) 

[ ~{n + l)(n + a + l)(n + 6 + l)(n + a + b + 1) 
X V (2n + a + b + l)(2n + a + b + 3) ' 

For a normalizable solution we must have (^1^) = 
J2^Lofn( e ) finite, thus f n should converge to zero. A 
closed-form expression for this sequence cannot be de- 
termined. One can observe though that for large n it 
asymptotically approaches 

/„ ~ C+ X* + C- A" with A± = , (9) 

C± being constants that depend on the minimal length 
through £ and the energy eigenvalues through e. This 
allows one to determine numerically the Coulomb spec- 
trum, by imposing C+ = 0. 

As an independent check, we used two different algo- 
rithms. First, for fixed minimal length £, we imposed 
fn+i/fn = A_ for sufficiently large n and scanned for 
the values of e for which the recursion (JSJ) gives /_i = 0. 
The contents following (J5J , and concluding with the first 
algorithm just described, represent results in unpublished 
work of Joseph Slawny. We thank him for making these 
available to us prior to publication. 

The second algorithm is more direct: for a given min- 
imal length £, we determined the values e for which /„ 
converges to 0. The subtlety is that C+ will never be 
represented internally as exactly zero, and the term cor- 
responding to it will eventually dominate our sequence. 
One can still identify the energy eigenvalues from the 
sign switch which occurs in C+ and correspondingly in 
the large-n behavior of /„. 

The algorithms yield consistent results, sampled in 
Fig. n We can see that the degeneracy among different 
angular momentum states is lifted: higher-^ states get 
smaller corrections. The only exceptions are the S states 
for j3' > 2(3 and D = 3. These states start out with 
the lowest, negative correction for small £, but cross the 
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FIG. 1: Plot of energy eigenvalues of the Coulomb potential 
in units of the regular ground-state energy as a function of the 
minimal length in units of the Bohr radius £ = hs/ (3 + /3'/oq. 
For principal quantum numbers n = 1, 2, 3, two extreme cases 
are represented: f3' = (thick lines) and j3 = (thin lines). 
Continuous line, t = 0; dotted line, t = 1; dash-dotted line, 
1 = 2. 

higher-angular-momentum levels as £ increases. Another 
important remark, expected but not readily transparent 
in the representation used, is that the energy values con- 
verge to the usual result in the limit when the minimal 
length is taken to zero. 

In order to get better insight into the observed behav- 
ior one can take another approach. This is also needed 
because for the very small values of the minimal length 
that interest us, i.e., several orders of magnitude below 
the Bohr radius, the numerical convergence becomes very 
slow and prone to rounding errors. 

As mentioned before, the choice of momentum repre- 
sentation, while convenient, is not necessary. In relation 
© one can use the "pseudoposition" representation 

h d 

Xi=Xi, Pi =-—— 10) 
i ox i 

in which the position operators Xi are not diagonal, ex- 
cept for the limit /3 = ft = 0. In this representation one 
can treat the f3- and /^'-dependent terms as small, and 
use perturbation theory to deduce the corrections to the 
energy spectrum. 

The operator R 2 can be written as 

R 2 = f 2 + e(R 2 )e+t; 4 (R 2 h^ (ii) 
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where the first-order correction acts on the radial part of 
the wave-function through 

+ „(,>_!)_ 3(D + l)]^L 

+ V (2L 2 + D 2 -D)-D{D + l), (12) 

with p = r/a and L 2 = £(£ + D — 1). 

In the expansion of the inverse distance R^ 1 = f^ 1 + 
£ 2 (-R _1 )£2 + 0(£ 4 ), w e expect on dimensional grounds 
that the first-order correction is a linear combination 
of terms of form (1/ p)d pp , (1/ p 2 )d p , and 1/p 3 . Indeed, 
substituting this form into Rr 1 R 2 R^ 1 = 1 determines 
uniquely 

(R- 1 )^ _ [(2D - 5) - T)(2D - 3)](D - 1) - A-qL 2 



+ 



4p 3 

(3-r))(D-l) d 1 d 2 



2p 2 



dp p dp 2 



(13) 



Thus, expressing the expectation value ((l/p)d pp ) 
through the Schrodinger equation and ((l/p 2 )<9 p ) in 
terms of (l/p 3 ), the first-order corrections to the energy 
eigenvalues can be written as 



e 



(Z>- 1) (377 -1 



+2 



where n 



1 /l 



n 2 \ p 



-i{i+i)(i- v ) 

(D-l)(l- V )p D - 3 



(14) 



D-3 



£ 



£ + and U n£ (p) is the 

unperturbed Coulomb radial wave-function. 

A note of caution is needed here. While the expansion 
is apparently in £ 2 , a quick calculation of higher-order 
terms confirms what is expected on dimensional grounds, 
namely, that the expansion parameter is /3/r 2 = £ 2 /p 2 . 
The £-quartic part in the expansion l|llf) of R 2 contains 
terms of the type £ 4 / p 2 and {£ 4 /p)d p . Therefore, the 
approximation 113|) for R^ 1 is no longer good for p < 
£. In particular, in the actual operator R^ 1 there is no 
singularity at the origin. 1 

Let us estimate the error. The largest discrepancy be- 
tween (|14|) and the actual value comes from the expec- 
tation value of 1/p 3 calculated over the interval [0,p c £] 
on which the approximation (|13|l breaks down, where 
p c = 7'c/ao ~ 1- F° r a n angular momentum state £, 
this is of the order 



t2l+D-\ 



2l+D-i 



dr 



for £ > or D > 3. (15) 



1 This is quite general. Any P-dependent commutation relation is 
expected to expand like and to exhibit this behavior. 
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FIG. 2: Comparison of different results for (a) 2P states with 
r\ = 0, 1 and (b), (c) IS states with r\ — 0, 1/3, 1. Solid lines, 
numeric result; dotted lines, perturbative result; dash-dotted 
lines, Ref. double-dashed lines, Ref. 0. The perturbative 
expression is given for (a) by 116H and for (b), (c) by formula 
111711 . with the coefficient C chosen such that it agrees with 
the numerics at £ 2 = 10 -6 . 



For D > 3 or £ ^ 0, this contributes only a higher- 
order term; thus it is safe to use (f 141) . and we finally 
arrive at 



En 



2? 



(£>-!) (377-1) 
U{£+l){£+\) 



n - 



l 



+ 



(16) 



This expression generalizes the result of Ref. [|| for 
arbitrary 77 and D. For the particular case D = 3 and 
77 = 1/3 (i.e., (3 1 = 2(3) it reduces to the one obtained 
there. Moreover, it is in excellent agreement with our 
numerical results. [See Fig. Ela) for the case of the 2P 
level.] 

When £ = and D = 3, the integral Ijl5|l is infinite. 
We can use only the part of (|13fl that can be trusted, i.e., 
we have to cut off the expectation value integral (l/p 3 ) 
at p = p c £. The leading-order terms are 

^ = ^^eEM) + ce + oi e), m 

ho n a 

where £d(p c £) = — ln£ — (7 + lnp c ) + 0(£) is the expo- 
nential integral function. The coefficient of the £ 2 term 
gets contributions from several parts. First, there are 
the remaining terms in l|14|) . Second, the actual value 



4 



of R^ 1 is bounded on the interval [0,p c £], so by cutting 
off the integral at p c £ we are neglecting another term of 
order £ 2 . Lastly, the exact choice of the cutoff value p c 
contributes another £ 2 term. Because we do not know 
the exact form of R , we cannot calculate analytically 
the second of these contributions. When needed, C can 
be determined numerically, by fitting relation (|17f) to the 
numerical results at a sufficiently low value of £ 2 . 2 

When compared to the numerical results, the behavior 
of the energy as a function of minimal length is nicely 
reproduced in Figs. Etb) and|2^c). Our results disagree 
with Ref. 0. The difference is well explained by the ne- 
glect there of all but linear terms in (3. These terms crit- 
ically affect the small- r behavior of RT 1 , and cannot be 
neglected. Reference Q arrives at a different expression, 
which is independent of r\. However, we could neither 
account for the discrepancy nor reproduce those results. 

We can finally set out to determine the constraint 
on the minimal length Aa; m j n from precision hydrogen 
spectroscopy. A naive estimate, obtained by imposing 
that the corrections are smaller than the experimental 
error on the value of the hydrogen 15-25 splitting, gives 
Ax min > 300 GeV (cf. H0). Unfortunately, this esti- 
mate would be correct only if the measured value of the 
physical observable agreed with the theoretical predic- 
tion and the main source of error were the experimental 
one. This is certainly not the case for the 15-25 split- 
ting in hydrogen: known to 1.8 parts in 10 14 , it is one 
of the most precisely measured quantities today and is 
considered a de facto standard |8( • The value of the Ry- 
dberg constant is determined using this measurement as 
an input, and thus the theoretical uncertainty is orders 
of magnitude above the experimental one. 

A better estimate is obtained by including contribu- 
tions of the (hypothetical) minimal length in the Lamb 
shifts. The strongest constraint is expected from the 15 
Lamb shift, being the one determined most precisely and 
getting the largest correction. The measured 15 state 
hydrogen Lamb shift of L°* pt = 8172.837(22) MHz @ is 
larger than today's best theoretical prediction L^g eor = 
8172.731(40) MHz [l(j by about ha experimental uncer- 
tainty. 

If we attribute the discrepancy entirely to the minimal 

2 For D > 3, (l/p 3 ) integrals are convergent, and the approximate 
solutions then agree nicely with the corresponding numerical re- 



length correction to the 15 state, the bound as a function 
of rj 1 obtained using the first two terms in (|17fl . is shown 
in Fig. [3J It is 1.75 GeV for r\ = 1/3 and increases to 
6.87 GeV for r\ = 1. Below r\ — 1/3, the constraint 
relaxes rapidly. Indeed, in this case the leading-order 
term in l|17fl is negative, and only the contribution from 
the next term can account for the observed difference. 
As a comparison, including only the leading term, we 
can obtain a bound only for rj > 1/3, with consistent 
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FIG. 3: Constraint on the minimal length obtained as a func- 
tion of -q, including the two highest-order terms (continuous 
line) and just the leading-order term (dash-dotted line) . 



results for r\ > 0.5. 

We should point out that the theoretical Lamb shift 
predictions are somewhat frail because of the uncertain- 
ties in the proton charge radius 1 1 11] . These are the same 
order of magnitude as the ones discussed here; thus one 
should consider the values in Fig.|21rather as upper limits 
for the minimal length. There is also the possibility of 
using muonium spectroscopy, but the current limits are 
still weak for our purposes. Details and other implica- 
tions for QED are under investigation. 
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